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Behavior of Charmonium systems after Deconfinement 
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We present a study of charmonia in hot gluonic plasma, for temperatures upto three times the 
deconfinement transition temperature Tc. qq systems with quark masses close to the charm mass 
and different spin-parity quantum numbers were studied on very fine isotropic lattices. The analysis 
of temporal correlators, and spectral functions constructed from them, shows that the J/^p and r]c 
survive up to quite high temperatures, with little observable change up to 1.5 Tc, and then gradually 
weakening and disappearing by 3 Tc. For the scalar and axial vector channels, serious modifications 
are induced by the hot medium already close to Tc, possibly dissociating the mesons by 1.1 Tc. 
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I. INTRODUCTION 

The behavior of strongly interacting matter in a hot 
and dense environment has been a subject of considerable 
theoretical and experimental research. At very high tem- 
peratures and densities, hadronic matter is expected to 
undergo a phase transition (or crossover) to a deconfined 
plasma state. Dedicated heavy ion collision experiments 
are aiming to create this state, the study of which will en- 
rich our understanding of strong interactions and which is 
also directly relevant for the early universe. The supres- 
sion of J/ip peak in the dilepton spectrum is one of the 
most important signals of this phase transition. Unlike 
light quarks, charmonia may survive as bound states even 
in a deconfined medium, due to Coulomb attraction be- 
tween the quarks. However, based on non-relativistic ar- 
guments, Matsui and Satz predicted that already at tem- 
peratures close to Tc, binding between quarks is reduced 
enough to dissolve J/ip, and proposed its suppression as 
a signal of the deconfinement transition P] . Several later 
studies, based on potential model calculations, predicted 
a pattern of dissolution, with the higher excitations dis- 
solving earlier, and J/ijj dissolving at a temperature ~ 
1.1 Hi. 

Since the validity of such potential model calculations 
at high temperatures are not a priori clear, a more re- 
liable way of studying properties of charmonia in a hot 
medium is clearly desirable. One possibility is to directly 
analyze the thermal Green functions of the corresponding 
states. For a long time it has been thought that numeri- 
cal study on the lattice is not a practical tool for calcu- 
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lations of dynamic properties of QCD at finite tempera- 
ture. However, it was shown recently that meson spec- 
tral functions, which are directly related to the real time 
correlators, can be extracted from temporal correlators 
calculated on the lattice using the Maximum Entropy 
Method (MEM) Some earlier investigations of me- 

son spectral functions (without the use of MEM) were 
presented in 0. The method has been successfully ap- 
plied at zero temperature 0, 0, S 111 ■ Additional difficul- 
ties are present at finite temperature because of the finite 
physical extent of the Euclidean time direction. However, 
the methods of Ref. were applied successfully to gain 
qualitative information about the medium modification 
of mesonic states with temperature 0,0,^3. Such stud- 
ies have provided useful quantitative information about 
dilepton rates j ^.nd quite unexpected nontrivial struc- 
ture of the low energy spectral function for the ss mesonic 
states Tl^. 

First applications of such methods to charmonium sys- 
tems also produced interesting and unexpected results 
|l2l Il3| |. The ground state charmonia, rjc and J/^, were 
found to survive well after the deconfinement transition, 
at least upto temperatures of 1.5 Tc [T^ . and no signifi- 
cant modifications of their masses were found on crossing 
the transition temperature Both these features are 
in sharp contrast to the existing potential model studies. 
The IP states Xc ^^nd xj, on the other hand, were found 
to undergo serious temperature modifications, possibly 
dissolution, already very close to Tc Later studies 

confirmed these results, and also found evidence of dis- 
appearance of the IS states from the spectral function at 
higher temperatures [T^[l5| . 

In the present paper we are going to report a de- 
tailed study of finite temperature charmonium correla- 
tors and spectral functions in quenched QCD using very 
fine isotropic lattices, expanding on preliminary results 
published in 0, ^|. The plan of the paper is as fol- 
lows. In the next section, we give a short discussion of 
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the definition of the spectral function and its properties. 
Sec. mil gives details of our lattices and simulation pa- 
rameters. In Sec. IIVI we explain our analysis methods, 
giving a short outline of the maximum entropy method. 
We then present the results below the deconfining tran- 
sition, which allows us to explain some features of our 
analysis that are used later. Sec. El is the central part of 
the paper, where we present our results for the charmo- 
nium correlators and spectral functions above Tc- Study 
of the systematics is a crucial part of the MEM analy- 
sis, at least at the current stage of the finite temperature 
mesonic spectral function studies. Therefore, a detailed 
discussion of the dependence of the results of Sec. El 
on different systematics follows in Sec. I VII In the next 
section, we discuss the spatial correlators and screening 
masses of charmonia, which show some interesting ther- 
mal effects and lead us to in-medium modifications of 
the dispersion relations for the J/V'- Finally, Sec. IVIIII 
contains a summary of our results, along with their phe- 
nomenological implications. We also shortly discuss here 
some points regarding the potential model studies. Read- 
ers who are interested in heavy ion phenomenology but 
not particularly interested in the details of our analysis 
may concentrate on Sections IVl and IVIIII 



II. SPECTRAL FUNCTION AND STATES AT 
FINITE TEMPERATURE 

Most dynamic properties of thermal systems are incor- 
porated in the spectral function. The spectral function 
o'h{po,p) for a given mesonic channel H in a system at 
temperature T can be defined through the Fourier trans- 
form of the real time two point functions and D^, or 
equivalently, as the imaginary part of the Fourier trans- 
formed retarded correlation function [la , 
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The correlators D^'^^\xo,x) satisfy the well-known 
Kubo-Martin-Schwinger (KMS) condition 



D>{xo.,x) = D<{xQ + i/T,x) 



(3) 



Inserting a complete set of states in Eq. ^ and using 
Eq. ©, one gets the expansion 
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A stable mesonic state contributes a 5 function-like 
peak to the spectral function: 
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where Z is the partition function, and fc"'™^ refers to the 
four-momenta of the state \n{m)). 



<yH{po,P) = mJH\H)\h{po)S{p'' - ml) 



(5) 



where mn is the mass of the state. For an unstable parti- 
cle one gets a smoother peak, with the peak width being 
related to the decay width. For sufficiently small decay 
width, a Breit-Wigner form is commonly used. As one in- 
creases the temperature, due to collision broadening the 
contribution of the states in the spectral function change, 
and at sufficiently high temperatures, the contribution 
from a state in the spectral function may be sufficiently 
weakened and broadened that it is not very meaningful 
to speak of it as a resonance any more. Such change 
of contributions of resonance states and eventual 'disap- 
pearance of resonances' in the thermal spectral function 
has been studied analytically, for example, in the Nambu- 
Jona-Lasinio model in Ref. \it\ . Finally, at very high 
temperatures one would expect the spectral function to 
consist only of a smooth continuum starting at twice the 
charm quark mass. The spectral function as defined in 
Eq. can be directly accessed by high energy heavy 
ion experiments. For example, the spectral function for 
the vector current is directly related to the differential 
thermal cross section for the production of dilepton pairs 
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The presence or absence of a bound state in the spectral 
function will manifest itself in the peak structure of the 
differential dilepton rate. 

In finite temperature lattice calculations, one calcu- 
lates Euclidean time propagators, usually projected to a 
given spatial momentum: 

Gh{t,p)^ I d\e'P-^{TrJH{r,x)JH{0,0))T (7) 

where (...)t means thermal average, and means or- 
dering in Euclidean time r. This quantity is the analytic 
continuation of D'' (xQ.p): 



GH{r,p) = D>{-tT,p}. 



(8) 



Using this equation and the KMS condition one can easily 
show that Gh{t, p) is related to the spectral function, Eq. 
by an integral equation: 



Gh{t,p) = 



dhjaH{uJ,p)K{uj,' 



cosh(w(T- l/2r)) 



(4) 



sinh(w/2r) 



(9) 
(10) 



Inverting Eq. @ one can extract spectral functions and 
properties of hadrons from correlators calculated in lat- 
tice QCD. In what follows, we use Eq. Q to extract the 
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behavior of degenerate heavy meson systems in a ther- 
mal medium from finite temperature mesonic correlators. 
Eq. H10|l is valid only in the continuum. It is not clear in 
general whether G{t,p) calculated on the lattice will sat- 
isfy the same spectral representation, but it was shown 
in Ref. that this is the case for the free theory. 



III. DETAILS OF THE LATTICE AND 
SIMULATION PARAMETERS 

In this work we restrict ourselves to the quenched ap- 
proximation and use pure gauge lattices generated using 
the isotropic Wilson action. In order to have enough 
points in the temporal direction at high temperatures, 
we need very fine lattices. We use lattices at three dif- 
ferent lattice spacings in the range ~ 0.02 — 0.05 fm. At 
the higher temperatures only the finer lattices are used, 
while at the lower temperatures, comparison between re- 
sults from our finer and coarser lattices gives an idea 
of the effect of the limited number of data points in our 
analysis. The lattices were generated using the heat-bath 
and over-relaxation algorithm, with each sweep consist- 
ing of one heat-bath step followed by 4 over-relaxation 
steps. The configurations were separated by 200 - 800 
such sweeps, the separation in each case being 5-8 times 
the autocorrelation time of the Polyakov loop. The ex- 
act details of the lattices are given in Table A subset 
of these lattices were used for the light meson studies of 
Ref. il^. 



TABLE I: Lattice parameters. The lattice spacing is obtained 
from the string tension. 





a-i(GeV) 


csw 




Size 


T/Tc #conf 


6.499 


4.04 


1.494176 


0.13558 


48=* 


X 24 


0.62 


50 










48=* 


X 16 


0.93 


50 










48=* 


X 12 


1.24 


50 










48=* 


X 10 


1.49 


45 


6.64 


4.86 


1.457898 


0.13495 


48=* 


X 24 


0.75 


100 










48=* 


X 16 


1.12 


50 










48=* 


X 12 


1.5 


60 


7.192 


9.72 


1.35500 


0.13437 


40=* 


X 40 


0.9 


85 










64=* 


X 24 


1.5 


80 










48=* 


X 16 


2.25 


100 










48=* 


X 12 


3.0 


90 



For the valence quarks, we use the clover-improved |23| 
Wilson action, taking the clover coefficients from the non- 
perturbative estimates of the ALPHA collaboration[2l] . 
The clover coefficients used, and the critical coupling Kc 
at each /3, are shown in Table |l| For our coarsest lattice, 
we use three k values that bracket the charm, giving a 
pseudoscalar mass in the range 1.7-4 GeV. This allows 
us to study the mass dependence of the properties that 
we are considering. At the finer lattices, we use one k 
value close to the charm mass. 



In order to study the lowest states in each of the four 
channels, we take operators with four different spin struc- 
ture. The operators we study are listed below, along with 
the names and spectroscopic representations of the lowest 
states for each channel: 



Jh 



cc 






C75C 


'So 


Vc 




'S^ 






'Pi 





(11) 



For zero temperature spectrum studies, smeared oper- 
ators are usually used to increase the overlap with the 
ground state, so one can study the ground state proper- 
ties already from correlators at small distance. However, 
when one wants to study other properties of a channel, 
for example meson decay constants, one has to be care- 
ful about using smeared operators. Here since we will be 
interested mostly in the existence versus dissolution of 
the states, it is tricky to use smeared operators. Smear- 
ing mimics bound states even when there are none j^, 
and it is also difficult to extract phenomenological quan- 
tities, like the dilepton rate in Eq. ©, if one works with 
smeared operators. Therefore we use point-to-point cor- 
relators in this study. When one knows that there are 
bound states, it may be better to use the smeared oper- 
ators to calculate the mass of the states. 

The use of the non-perturbative clover action removes 
0(a) discretization errors. However in the case of heavy 
quarks discretization errors of order 0{am) can be large. 
The errors in using the clover quarks due to finite mass 
of the quarks have been discussed in detail in Ref. [2^ . 
The main sources of 0{ma) errors are in the difference 
between the pole mass and the kinetic mass, and in the 
renormalization constant used to connect the lattice op- 
erators to continuum ones. For finite ma, the pole mass 
that controls the fall-off of the correlator differs by 0{ma) 
terms from the physically important kinetic mass that 
controls the dispersion relation. This can already be seen 
from the free quark dispersion relations |22j |: 



amo 
ampoic 



1 1 

log(l + amo) 
amo(l + amo)(2 -t- am,o) 
2 -I- 4amo -I- mn 



(12) 



where mp is the standard definition of the bare mass, 
Wpoie governs the fall-off of the free quark correlator 
and TOkinotic is the term governing the quadratic momen- 
tum dependence in the dispersion relation, and therefore, 
controlling much of the spectrum and other interesting 
physics in heavy quark systems. A large mismatch be- 
tween mpoio and mkinotic indicates that the quark is too 
heavy to be treated relativistically. Our lattices are fine 
enough that in most cases the error due to the finite mass 
of the charm quark is small enough. The hopping param- 
eters we have used for our different lattices are listed in 
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Table [HJ We also list there the tree level pole and ki- 
netic masses, using Eq. H12|l . The tree level mismatch 
between mpoic and mkinetic is < 5% for all the sets, ex- 
cept the heaviest quark for set I. To see whether the error 
is also small for the interacting theory, we also looked at 
the pole and kinetic masses of the mesons, where the pole 
mass is obtained from the zero momentum correlator and 
the kinetic mass is obtained from the dispersion relation 

E'^Ml,^ + ^p\ (13) 

While Mkin is quite noisy, we found that there is hardly 
any significant deviation from the relativistic dispersion 
relation for the quark masses used by us, except possibly 
for set IC ^3 ■ We therefore use a relativistic treatment 
for our quarks in what follows, and set the masses of the 
mesons from the pole masses. 

TABLE II: Hopping parameters used and the pole and kinetic 
masses for different sets. 



13 




arripoie 


ffl'Tlkinctic 


amo 


Tag 




0.1325 


0.0889 


0.0893 


0.0860 


I A 


6.499 


0.1300 


0.1584 


0.1608 


0.1583 


I B 




0.1234 


0.3328 


0.3531 


0.3640 


I C 


6.64 


0.1290 


0.1800 


0.1835 


0.1821 


II 


7.192 


0.13114 


0.0940 


0.0945 


0.0916 


III 



The lattice operators are connected to the continuum 
operators as 

Jg™* = 2KZH{a, m, ^ = 1/a) jjj^V-^ (14) 

The renormalization factors ZHicL,m, fi — 1/a) are esti- 
mated using 1-loop tadpole improved perturbation the- 
ory 

Znia, m, fi = 1/a) = uo{a){l+av{l/a)zH){l+bH{a)ma) 

(15) 

bnia) = (1 + av{l/a)B„)/un. (16) 

Here ay (1/a) is the running coupling calculated from 
the heavy quark potential at scale jj. — 1/a (2J]. The 
one loop coefficients zh and bn were calculated in Refs. 
psj l and [2^ respectively. For the scalar and pseudoscalar 
channels zh was evaluated at scale ji = 1/a. We evaluate 
Wo from the measured plaquette, 

Mo = (gTrf/plaq)^. 

The resulting renormalization factors are shown in Table 

inii 

We close this section with a discussion of the zero tem- 
perature masses of the charmonia corresponding to the 
K values given in tabled For our coarsest lattices, we 
use three k values that bracket the charm quark. For the 
other lattices, we use one k, value each. Since we have not 
performed calculations on symmetric (zero temperature) 



TABLE III: Renormalization factors of local operators for dif- 
ferent sets. 



p 




Zps 


Zsc 


Zvc 


Zax 




0.1325 


0.782 


0.839 


0.900 


0.926 


6.499 


0.1300 


0.847 0.913 


0.975 


1.003 




0.1234 


1.032 


1.124 


1.188 


1.221 


6.64 


0.1290 


0.881 


0.948 


1.007 


1.034 


7.192 


0.13114 


0.839 


0.886 


0.936 


0.957 



lattices, our estimates of the zero temperature masses 
are based on exponential fits to spatial correlators. In 
fact, if one considers the spatial correlators on lattices 
below deconfinement, the finite temperature lattices can 
be considered as small (in one of the directions) zero tem- 
perature lattices. For the lattices below Tc used in this 
study (see TablePi, the small extent is about 1 fm for the 
two coarser lattices (sets I and II) and about 0.8 fm for 
the finest lattice (set III). The effect of finite volume on 
the masses of charmonia was studied in detail in j23| and 
it was found that even for lattices of size 0.75 fm there 
are no sizeable finite size effects. Furthermore, even for 
the much larger light mesons, spatial correlators calcu- 
lated on lattices at 0.9 Tc have been seen to give good 
estimates of zero temperature masses The masses 

we get from exponential fits are given in Table HVl 

TABLE IV: Masses (in GeV) of the different charmonium 
states obtained from one exponential fits of spatial correlators 
below Tc- 

Set Mass (in GeV) 

»?c J/ip XcO Xcl 

I A 1.673(8) 1.806(8) 2.151(27) 2.299(44) 

I B 2.440(5) 2.517(10) 2.874(30) 3.006(31) 

I C 4.151(4) 4.191(6) 4.600(44) 4.670(48) 

II 3.071(6) 3.141(12) 3.526(53) 3.622(117) 

III 3.744(14) 3.803(10) 4.442(86) 4.666(174) 

For set II, which is very close to the physical value of 
the charm quark mass, we get 70 ± 12 MeV for the hy- 
perfine splitting, consistent with the analysis of Ref. |27l |. 
The mass splitting between P and S states is also consis- 
tent with the findings of j23|. For our finest lattices the 
splitting between P and S states is largely overestimated. 
This could be due to the limited (~ 0.82 fm) extent of 
the lattice. 



IV. MESON CORRELATORS AND SPECTRAL 
FUNCTIONS FOR T < 

We begin our discussion of spectral functions and tem- 
perature effects on meson properties with an analysis of 
the meson correlators below Tc. We look at the temporal 
correlators, Eq. Q, for the four channels in Eq. \11\ 
and for the lattices with T < Tc as specified in Tabled 
For the lattices at /3 = 6.64 and 7.19, we have one set 
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each below Tc, the spatial correlators from which were 
used in the previous section for setting the parameters. 
Here we use the temporal correlators and investigate the 
spectral function at temperatures T <Tc. On our coars- 
est set dX (5 — 6.49, we have lattices corresponding to two 
different temperatures below Tc- This allows us to check 
for any change in meson properties as one approaches Tc 
from below. For the most part of this study we use zero 
momentum projected correlators only, and refer to them 
as G(t), and the corresponding spectral function as <j{uj) 
(we also suppress the subscript H from now on). 

The operators corresponding to the different chan- 
nels have been listed in Eq. CIJ- For the vec- 
tor channel, we take the trace over all four directions, 
S^=o ^^'^ *° current conservation, the time 

component, Gqq^{t), is a constant, and therefore, con- 
tributes only a delta- function at ld — Q to the spectral 
function. In the free quark case, this delta function can- 
cels a similar contribution from the space-averaged cor- 
relator [3ll | . In the axial vector channel, we sum over the 
three spatial directions, X]i=i Gf^^ [t). 

At least in principle, one can obtain all the informa- 
tion about the spectrum in a given channel by inverting 
Eq. lO to calculate (j{lo). Since our lattices have only 
0(10) independent data points in the temporal direction, 
however, extracting a{uj) from lattice data for G(t) is 
a highly nontrivial problem. Considerable progress to- 
wards the solution of this problem has been made in the 
last couple of years by using Bayesian methods . In the 
next subsection, we explain the rudiments of the maxi- 
mum entropy method (MEM), as employed in our anal- 
ysis; in the rest of this section, we use this method to 
extract the spectral function, and discuss the results. 

A. Analysis; MEM methods 

Bayesian techniques for extracting information from 
inadequate data depend on providing prior information 
to the analysis in some form j^^. In the maximum en- 
tropy method (MEM) of data analysis, the prior is usu- 
ally introduced in the form of an 'entropy term' by using 
the positivity of the spectral function (g(a;) > 0,w > 0) 
to give it a probability interpretation |29| . Quite general 
principles lead to the Shannon entropy 

/>oo 

S= duj[a{Lu) — m{Lu) ~ a{uj)log{a{uj)/m{uj))] (17) 

^0 

(see Refs. |^ for justifications for the form of the 
entropy term). The entropy depends on an arbitrary 
function m{uj), which is called the default model. This 
function can be considered as a part of our prior infor- 
mation. The central problem becomes then to find the 
a{Lu) that maximizes the 'free energy' 

F = L-aS (18) 

or equivalently the conditional probability P[a\DH] of 
having the spectral function a given some data D and 



the prior knowledge H |2^. Here L is the likelihood 
function for the data D, i.e., L — and a is a real 

and positive parameter. 

At zero temperatures, MEM has been successfully ap- 
plied in spectrum calculations 0, 0, Q, 13 ■ At finite tem- 
peratures, where one is limited further by the small phys- 
ical extent of the temporal direction, precise quantitative 
information is more difficult to obtain. Still, useful qual- 
itative information has been obtained about the shape of 
the spectrum and the thermal dilepton rate 'ol, 0, ^| . 

We follow here Bryan's algorithm 30] for inverting Eq. 
@. This method uses a singular value decomposition of 
the kernel, to find the relevant, nonsingular directions in 
the inversion process, which are less than the number of 
data points. Then the search for solutions is restricted 
to this space, which is now a well-defined problem. Also, 
in this algorithm one avoids any explicit dependence on 
the external parameter a in Eq. (|18|) by integrating over 
a with a conditional probability (see ^, ^] for details). 

To reconstruct the spectral function the default model 
171(0}) needs to be specified. Since we want to explore the 
modification of the ground state peaks in the spectral 
function, it is natural to choose the default model to de- 
scribe only the high energy part of the spectral function. 
Due to asymptotic freedom the spectral function at very 
high energies is expected to approach the case of two free 
quarks. For two free quarks in a mesonic channel H, the 
spectral function in the continuum is given by 31] 

aniuj) = —LU^Qiuj^ - Am^)t(inh(uj/4T)Jl - 4m^/uj^ 

[an + —bn] (19) 

3 

« — jflHi^^ for ui — > 00. (20) 

OTT 

For the scalar, pseudoscalar, vector and axial vector 
channels, the constants an are 1, 1, 2 and 2, respec- 
tively 31]. On the lattice, the spectral functions for the 
different channels have been calculated recently for the 
free theory ^9,] . They show a considerably different high 
energy structure from Eq. H20|) . Interaction introduces 
additional modification of the large oj behavior of the 
spectral function j^. At low temperatures, where the 
temporal extent is relatively large and the dependence on 
the default model is less important, we use m(u;) = mow^ 
for the default model, with mo — 3/(87r^) for scalar and 
pseudoscalar channels and toq — 3/(47r^) for vector and 
axial vector channels. Lattice introduces an ultraviolet 
energy cutoff 

2 

Wmax = - ln(7 -I- amo) w 4a 

for our masses. This is the frequency cutoff we use in our 
analysis. Above Tc the use of a more refined form of the 
default model, with more accurate information about the 
high energy structure of the interacting theory, becomes 
necessary, which we discuss in the next subsection. 
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Use of the full covariance matrix is desirable in the 
MEM analysis. However, it has been shown that with 
limited statistics sometimes the small eigenvalues of the 
covariance matrix can be too noisy and it may be better 
to use some smoothening of the smaller eigenvalues [33| . 
In our analysis, we use the full covariance matrix in most 
cases; but sometimes the use of the full matrix gives rise 
to an unphysical peak at w ^ 0. In such cases, we use 
the minimum smoothing that can get rid of this peak, if 
it is possible with a small amount of smoothing. 

For the spectral density extracted from noisy data with 
MEM, one can assign error bars to the spectral density 
integrated over small regions of frequency (see Ref. [2^): 
one defines an 'integrated spectral density', 

H{uj,A) = / duj(T{uj) 

and then the covariance of this quantity is calculated as 

{SH^) = J dcj J duj'a{Lo)(j{uj'). (21) 

Since we are concerned with the peak structures in the 
spectral function, in the following sections we show the 
error bars on a region integrated over the peak width 
of each peak, to show the statistical significance of the 
peaks. 

B. Discussion of the spectral functions 

As mentioned in the previous section, the reliability 
of the spectral function reconstructed from the temporal 
correlators using MEM depends on the number of data 
points. We start the reconstruction of spectral functions 
below Tc on our finest lattice: set III. Besides allowing 
us to have the maximum resolution, fine-graining of the 
lattice also has the advantage that the lattice artifacts 
are shifted further away from the ground state peaks, 
allowing an accurate extraction of the ground state peaks 
of interest. 

Fig. n shows the spectral functions at 0.9 Tc, recon- 
structed from the correlators for this set in the different 
channels. In order to remove the uninteresting uP' rise 
in the spectral function, here and in everywhere else we 
plot the 'reduced spectral function' 

p{u}) = a{u})luj^. (22) 

At r = 0.9Tc the physical distance available for the 
analysis of the temporal correlators is about 0.82 fm, and 
we have 40 data points. This allows us to have a good 
resolution of the ground state spectral function. We also 
found that omitting the data point at smallest r enhances 
the peak structure for this set. In each channel we find 
a three peak structure, with the ground state peak repli- 
cating quite well the properties of the zero temperature 
ground state in Table Hvl We also find that for energies 
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FIG. 1: Spectral functions at 0.9 Tc for the different channels, 
for set III. 

uj ~ 2/a, where the doubler states are expected to con- 
tribute, the structure of the spectral function is consid- 
erably different from the free spectral function. Similar 
peak like structures have been seen in Q . As we explain 
below, we expect the additional peak seen at w 1/a is 
also a lattice artifact. 

A comparison with the spectral structure obtained 
from set II helps in clarifying the nature of the second 
peak in Fig. ^ The spectral functions in different chan- 
nels for this set are presented in Fig. [51 For this figure we 
use again the free lattice spectral function as the default 
model. The properties of the zero-temperature peak in 
Table Hvl are reproduced quite well. The other peaks are 
seen to shift at approximately the same ratio from Fig. 
Has the inverse lattice spacing. This confirms that also 
the second peak obtained in Figs. ^IHis dominated by 
lattice artifacts. 
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FIG. 2: Spectral functions at 0.75 Tc for the different chan- 
nels, for set II. 

In our coarsest set of lattices, at /3 = 6.499, we have 
two temperatures below Tc- This allows us to study the 
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temperature dependence of the mesons below Tc, and will 
help in introducing certain elements in our analysis, as 
we will see below. 

For the correlators at 0.6 Tc the results are shown for 
set IB in Fig. |2| We get a three peak structure for 
the different channels, similar to the one shown in the 
previous figures. The ground state can be resolved quite 
reliably, and the strength and position of the ground state 
for the different quark masses agree with the ground state 
measurements obtained in previous section. 
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high Lu side is considerably broader. However, we expect 
this difference in the spectral function to be an artifact 
of the limited number of data points and physical dis- 
tance at this temperature. In fact. Fig. 21 does indicate 
that at 0.9 Tc, the spectral analysis can not resolve the 
three peak structure, and the ground state peak is con- 
taminated by the peak at w ~ 1/a. One way to improve 
the reconstruction of the ground state when the physical 
distance is small is to provide more accurate information 
for the high energy part in the default model. Since the 
physical distance covered is small, the high energy part 
is still important in the correlator and it is only when we 
provide suitable information about the high energy part 
that it is possible to extract the information about the 
ground state accurately. From the analysis of the spec- 
tral function at the lowest temperatures we see that the 
high energy part of the spectral function consists of two 
broad peaks. This is a general feature of the spectral 
functions which does not depend on the default model 
m{uj) as we will see in section IVll So we use a default 
model where the high energy part of the spectral function 
is taken from the spectral function constructed at 0.6 Tc, 
and in the low energy part we use the model 

ni{uj) — miuP' 



FIG. 3: Spectral functions for different channels at 0.6Tc for 
set IB. 



Free default, 0.61^ 
Free default, 0.9T. 
2-Pk default, 0.6T. 
2-Pk default, 0.9T. 




FIG. 4: Spectral functions for the vector channel below Tc, 
using both the free spectral function and the high energy part 
of the spectral function from interacting theory as default 
model. Set IB. 

Figure 01 compares the spectral functions for the vector 
channel at 0.6 and 0.9 Tc, obtained with the free spectral 
function as the default model. The high uj part of the 
spectral function at 0.9 Tc is seen to be considerably dif- 
ferent from that at 0.6 Tc- Also, while p{u}) shows a clear 
ground state peak, the peak position is slightly changed 
from that at 0.6 Tc and the fall-off of the peak in the 



with 7TV\ chosen such that the default model is contin- 
uous [13 • The spectral function constructed at 0.9 Tc 
using this default model is also shown in Fig. ^ where 
for comparison we also show the spectral function at 0.6 
Tc using the same default model. As we can see, at 0.6 
Tc the spectral function correctly produces the ground 
state peak, in both the position and the strength of the 
peak. (In our present analysis with limited statistics the 
width of the peak is not of direct physical relevance. It 
is possibly an effect of the finite statistics.) Using this 
default model, now, the ground state at 0.9 Tc is recon- 
structed accurately, with the position and strength of the 
peak similar to that at 0.6 Tc. This also shows the use- 
fulness of using such a default model, and we will use it 
for reconstructions of the spectral function above Tc in 
the next section. In essence, what is being done is to 
realize that the spectral function for very large uj » T 
is dominated by lattice artifacts, and does not change 
with temperature. To reconstruct the ground state accu- 
rately from correlators at short distances it is important 
to use correct information of this high energy part. An- 
other way to confirm that the spectral function does not 
show any significant change between 0.9 Tc and 0.6 Tc 
is by comparing the correlators at 0.9 Tc with the corre- 
lators reconstructed from the spectral function at 0.6 Tc 
Since the temperature dependence of the correla- 
tors in Eq. Q comes from both the known temperature 
dependence of the kernel, K{uj,t), in Eq. ifTUI) . and the 
nontrivial temperature dependence of the spectral func- 
tion, a(w,p), we can try to focus on the temperature 
dependence of the spectral function by comparing the 
correlator with the reconstructed correlator at that tem- 
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perature, 

Grccon,T. (t, T) ^ J dua{uj, T*)K{lu, t, T) (23) 

where a{Lu,T*) is the spectral function at another tem- 
perature T*. The comparison of the correlators for the 
different channels at 0.9 Tc with the correlators recon- 
structed from cr(w,0.6Tc) for the corresponding channels 
is shown in Fig. This figure clearly demonstrates that 
the correlators at 0.9 Tc are completely described by the 
spectral function at the lower temperature. 



function constructed at the lowest available temperature 
to reconstruct the correlator at higher temperatures, as a 
way of studying the temperature dependence of the spec- 
tral function. This often provides a more robust analysis, 
since the lowest available temperature provides us with 
the maximum physical distance, allowing for a more reli- 
able extraction of the spectral function. Both these tech- 
niques will be particularly relevant in the next section 
for studying the spectral function above Tc, where small 
physical distance makes the extraction of the spectral 
function even more difficult. 
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FIG. 5: Ratio of the measured correlators at 0.9 Tc with 
the reconstructed correlators at this temperature, using the 
spectral function at 0.6 Tc (see text), for different channels 
and set IB. For visual clarity, the vector and scalar channels 
have been shifted horizontally. 

The results from the other quark masses for set I are 
very similar and show that for the quark mass range ex- 
plored by us, the properties of the ground state mesons 
do not change at least up to a temperature of 0.9 Tc. 
This is consistent with other earlier quenched studies in 
the mesonic sector 28j and the gluonic sector , which 
showed that the properties of the mesonic and gluonic 
states do not change from their vacuum properties till 
quite close to the deconfinement transition temperature. 
Of course, this result is not unexpected since one does not 
expect substantial thermal excitation of the heavy glue- 
balls already at these temperatures. Presumably, such a 
result may not hold in full QCD, where excitation of pi- 
onic states may produce observable effects already away 
from Tc. 

To summarize, in this section we have studied the spec- 
tral functions below Tc for heavy degenerate mesons in 
different channels. The ground state properties can be 
suitably reproduced by the MEM analysis of the tem- 
poral correlators. We also show that using nontrivial 
realistic information about the high energy structure of 
the spectral function in the interacting theory allows a 
more accurate reconstruction of the ground state already 
from short distance correlators available on a small num- 
ber of data points. We also propose using the spectral 



V. CORRELATORS AND SPECTRAL 
FUNCTIONS FOR T > Tc 

In this section we present our results for meson cor- 
relators above Tc. We first present results for set II, so 
that we can start with results close to Tc, and then go to 
higher temperatures with set III. 

A number of interesting and important statements 
about system modifications of charmonium spectral func- 
tions can be made by inspecting the finite temperature 
meson correlators. As was discussed in the previous sec- 
tion, the change in a{uj) with temperature can be studied 
already from the correlators, by first taking out the trivial 
temperature dependence of the kernel if (w, r, T) accord- 
ing to Eq. We construct Grocon(T, T) using a{u!, T*) 
from the spectral functions obtained for this set at 0.75 
Tc (see Fig. The deviations of Grocon(T, T) from the 
correlators directly measured above Tc then indicate sys- 
tem modifications of the spectral function above Tc- This 
method is more robust as it avoids using MEM at the 
higher temperatures, where one has less temporal extent 
and less data points. 

In Fig. Ela) we show the results of G(r, T)/Grccon(T, T) 
for set II for the pseudoscalar and vector channels (the 
IS states). For the pseudoscalar channel, the ratio is 
seen to remain equal to one upto 1.5 Tc, indicating no 
significant system modification of its properties upto this 
temperature. This is already in sharp contrast to earlier 
potential model results, which predicted a dissolution of 
rjc at « l.lTc. For the vector channel, one finds that the 
reconstructed correlator explains the data quite well at 
small distances. Some significant system modifications, 
however, are manifested in the correlators at distances 
> 0.25 fm. At 1.5 Tc, some modifications are already 
seen at distances > 0.2 fm. Later in this section, we 
will discuss possible system modifications that can cause 
such a change. 

The correlators for the scalar and axial vector channels 
(the IP states), are shown in Fig. Elb). The situation is 
seen to be quite different here: already at 1.1 Tc, a signif- 
icant system modification of the mesons is manifestated 
in G(t, T)/Grecon(T, T) at all distances. The measured 
correlators are seen to be considerably larger than the 
reconstructed correlators. Such a behavior can be quali- 
tatively understood if the IP states are dissolved already 
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Tc, smoothly matched to toicj^ in the low energy region 
(see Sec. Wm . 



FIG. 6: Ratio of the measured temporal correlators with the 
correlators reconstructed from the spectral function at 0.75 Tc 
for Set II. a) Pseudoscalar (PS) and vector (VC) channels; b) 
Scalar (SC) and axial vector (AX) channels. For visual clarity, 
the PS and AX channel results have been shifted vertically. 



at these temperatures: at low temperatures the spectral 
function for these states is zero below ~ 3.5 GeV. If the 
state is dissolved above Tc, one would expect the cor- 
relator to pick up contributions at significantly smaller 
uj ~ 2mc, making it larger. 

The comparison with reconstructed correlators can 
only tell us whether any medium modification can be 
expected at a given temperature. The ratios of the 
correlators give a rough idea of the magnitude of such 
modification. To further explore the nature of the finite 
temperature spectral functions it is necessary to recon- 
struct the spectral function from the correlators, using 
MEM. Above Tc, the physical distance available to us 
is small, and fewer data points are available, making it 
difficult to reliably extract the spectral function without 
any prior knowledge. However, as we have discussed in 
Sec. IIVBI introducing information about the high en- 
ergy part of the spectral function into the MEM analysis 
in the form of an appropriate default model we can over- 
come this problem. Therefore, in the default model we 
use the high energy part of the spectral function at 0.75 
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FIG. 7: The spectral function for the a) scalar and b) axial 
vector channels constructed from the temporal correlators, for 
Set II at temperatures between 1.12 and 1.5 Tc. The default 
model uses the high energy part of the spectral function in 
Fig. |5| as explained in the text. 

The spectral functions for the scalar and axial vector 
channels above Tc, reconstructed this way, are shown in 
Fig. [7| The significant ground state peak of Fig. [3 for 
these channels are not found at these temperatures. A 
non-zero spectral function is seen at a significantly lower 
UJ, but the peak structure is statistically not significant, 
and it could be related to a branch cut coming from two 
on-shell propagating quarks. Such a a{uj) at low uj can 
also explain the rise of the correlators seen in Fig. 
b). Since the quark mass for set II is very close to the 
physical charm quark (see Table lIVp , Fig. [7| together 
with Fig. Elb), will indicate that the Xc ^^'^ the xl are 
seriously modified, possibly dissolved, already at 1.1 T^. 

The situation is different for the IS states: FigureEla) 
seems to rule out dissolution or a major modification of 
these states. The reconstructed spectral functions, shown 
in Fig. |S1 support this: a statistically significant ground 
state peak is seen upto 1.5 Tc, showing that the rjc and 
the J/tp survive at least upto this temperature. While no 
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FIG. 8: The spectral function for the a) pseudoscalar and b) 
vector channels constructed from the temporal correlators, for 
Set II at temperatures between 1.12 and 1.5 Tc- The default 
model uses the high energy part of the spectral function in 
Fig. 13 as explained in the text. 



major change in mass is indicated, a reduction in peak 
strength is seen at 1.5 T^] however, as discussed below 
and also Sec. IVIL Fig. ^1 we expect it to be related to 
the associated systematics. 

Having seen that the J/ip survives till 1.5 Tc, it will 
be interesting to explore its properties at higher temper- 
atures. Our finest lattices at set III allow us to go upto 
a temperature of 3 Tc- The finer grid at 1.5 Tc here also 
allows one to check the reliability of the results of set 
II. FigureEla) presents G(t, T)/Grecon(T, T) for different 
temperatures. Grccon(''', T) is constructed using the spec- 
tral functions at 0.9 Tc presented in Fig. ^ Upto 1.5 Tc, 
the features of Fig. I^a) are very similar to that of Fig. El 
a). For the pseudoscalar channel, the data is completely 
explained by the spectral function at 0.9 Tc. Some small, 
but statistically significant, deviations appear at 2.25 Tc. 
At 3 Tc, we see significant modification at all distances. 
For the vector channel, at 1.5 Tc the data agree with 
the reconstructed correlator at small distances, but start 
showing significant modifications at larger distances, r > 
0.15 fm. At 2.25 Tc such modifications appear already at 



- 0.1 fm while at 3 Tc, G(r, T) /G,ocon(r, T) differ signif- 
icantly from 1 at all distances. 

For the scalar and axial vector channels, one sees at 1.5 
Tc a very similar enhancement of the correlator to Fig. 
Elb), with G(t, 1.5Tc) ~ 60% larger than Grocon(r, T) at 
^ 0.25 fm for the scalar and ~ 120% larger for the axial 
vector. The pattern of enhancement continues as one 
goes to higher temperatures. 
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FIG. 9: Ratio of the measured temporal correlators with the 
correlators reconstructed from the spectral function at 0.9 Tc 
for Set III. a) Pseudoscalar (PS) and vector (VC) channels; b) 
Scalar (SC) and axial vector (AX) channels. For visual clarity, 
the PS and AX channel results have been shifted vertically. 

To extract the spectral function from the temporal cor- 
relators above Tc, we use, as input to the default model, 
the large u structure shown in Fig. ^ Fig. ^ja) shows 
the spectral functions above Tc obtained with this default 
model. A strong and statistically significant ground state 
peak is obtained at 1.5 Tc. The peak position and peak 
strength are found to be very similar to those at 0.9 Tc, 
indicating essentially no significant change for the pseu- 
doscalar upto this temperature. This is completely con- 
sistent with the trend seen for G(t, T)/Grecon(''", T) for 
the pseudoscalar. At 2.25 Tc, we still see a statistically 
significant peak, but with a reduced strength, and at 3 
Tc the structure at the peak position is very weak and 
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not statistically significant. 
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FIG. 10: The spectral function for a) pseudoscalar and b) 
vector channels constructed from the temporal correlators, 
for Set III at temperatures upto 3 Tc- The default model uses 
the high energy part of the spectral function in Fig. Q as 
explained in the text. 

For the vector channel, the temperature dependence of 
the spectral function is very similar. A very strong peak 
is seen at 1.5 Tc, with no observed reduction in the peak 
strength from that below Tc. Unlike the pseudoscalar 
channel, we see here a slight shift of the peak position, 
but the amount of the shift depends on the systemat- 
ics and we are unable to make a definite comment on 
it. At 2.25 Tc, a significant peak is seen with a reduced 
strength, and at 3 Tc, the structure at the peak position 
is too broad to be interpreted as a resonance, and is also 
statistically not significant. 

As we will discuss in detail in the next section, the 
details of the features presented in Fig. [TUl and Fig. |S1 
are not dependent on the systematics incorporated in ex- 
tracting the spectral function, either from the specifics of 
the default model or from the small temporal extent and 
limited number of data points at higher temperatures. 
The IS state mesons survive with no significant change 
of strength at least upto 1.5 Tc, and at 2.25 a signifi- 
cant peak still survives but with a reduced strength, while 



at 3 Tc system effects reduce the peak strength enough 
that one should consider the state to be dissolved. Such 
a change in the spectral function is also completely con- 
sistent with the behavior of G{t,T) /Giccon{T,T) for the 
pseudoscalar channel, since a reduced ground state peak 
will cause a depletion of G(r). For the vector channel, 
the general picture indicated by Fig. ^|b) is the same. 
However, we also see here a depletion of G(r) at large 
distances already at 1.5 Tc, but no reduction of the peak 
strength in Fig. E|b), and a small increase in mass. 
While such an increase in mass at 1.5 Tc is consistent 
with the depletion of the correlator, we cannot isolate 
any physical mass shift from the systematics. 

We would like to comment here on similar works on 
J/ip in the quark-gluon plasma. Ref. |l3 uses smeared 
operators, which allows them to study the properties of 
the in-medium rjc and J/ip in more detail, and finds a 
large Breit-Wigner width 120 MeV and 210 MeV, re- 
spectively, at 1.1 Tc. They also find no reduction of the 
masses of the IS states above Tc. With our point-point 
correlators we cannot extract a width reliably; however, 
any such width will cause an enhancement of the cor- 
relator, which we do not see for these states. However, 
it is possible that such an enhancement is shielded by a 
corresponding reduction in peak strength. Ref. [14| uses 
point-point correlator, and concludes that the IS states 
dissociate already at 1.9 Tc- Of course, our quarks for set 
HI are a little heavier than physical charm, and that can 
cause the mismatch (the mass dependence of the mod- 
ification pattern is discussed below); however, we think 
another likely reason could be that due to the use of the 
free continuum spectral function, Eq. H19I) . as the default 
model in Ref. jl4j . (which does not provide accurate in- 
formation about the structure of the high energy part of 
the spectral function), and the small temporal extent at 
this temperature, MEM cannot resolve the peak struc- 
ture correctly. 

For our coarsest lattices, set I, we have three quark 
masses, giving pseudoscalar masses in the range 1.7 to 
4.1 GeV (see Table HVj). This allows us to study the 
mass dependence of the dissolution pattern, i.e., to an- 
swer the question whether the "melting temperature" for 
the quarkonia has a mass dependence. 

In Fig. ^Jwe show the ratio of the pseudoscalar corre- 
lators at different temperatures with the correlators re- 
constructed from spectral function at 0.62 Tc using Eq. 

At 0.93 Tc, expectedly, the correlators are explained 
by the spectral function at 0.62 Tc for all quark masses. 
For the lightest quark studied by us (k = 0.1325) some 
system modification is already seen at 1.24 Tc, and the 
modification becomes larger as one goes to 1.5 Tc. For the 
next heavier quark (k=0.13), at 1.24 Tc the correlators 
agree completely with the reconstruted correlators, and 
one sees temperature modification of the pseudoscalar 
only at the higher temperature. For our heaviest quark 
(k=0.1234). Fig. ^2 a) shows that the pseudoscalar cor- 
relator is completely explained by the spectral function 
upto at least 1.5 Tc. The vector channel shows a simi- 
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lar trend, with the lighter quarks showing larger medium 
modification at a given temperature. For the scalar and 
axial vector channels, we found significant medium modi- 
fications similar to Figs. Elb) and|51b) for all three quark 
masses, indicating that also for the heaviest quarkonia 
studied by us, the IP states undergo significant medium 
modification, possibly dissolution, already at 1.24 Tc- 
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FIG. 11: The mass dependence of G(r, T) /Grccon (r, T) at dif- 
ferent temperatures, for the pseudoscalar channel. The refer- 
nce spectral function used for 0^0011(1", T) is (j{uo, 0.62rc). For 
clarity, the points at 1.24 Tc have been shifted vertically, as 
indicated in the figure. 



model, with those shown in section IV where the mass- 
less continuum spectral function was used as a default 
model. For the quark mass in the free lattice model we 
use ma ~ 0.17, which corresponds to the bare mass for 
this set. As one can see from the figure the dependence 
on the default model is quite small, even though the two 
default models have very different structures. For com- 
parison, the masses extracted from the spatial correla- 
tors at this temperature are also shown as vertical lines. 
There is a small deviation between the peak position and 
the screening masses, ~ 3% for the pseudoscalar channel 
and ~ 7% for the vector channel. Since we do not ex- 
pect medium modification of quarkonia properties below 
deconfinement (cf. Fig. O this should be interpreted 
as systematic error in our MEM analysis. Similar state- 
ments also hold for the spectral functions at 0.9Tc. 




VI. DISCUSSION OF POSSIBLE SYSTEMATIC 
ERRORS IN MEM ANALYSIS 

The extraction of the continuous spectral function 
from 0(10) data points is, to start with, a rather del- 
icate problem, and one needs to understand the effect of 
the systematics well before drawing physics conclusions. 
In this section we are going to discuss possible systematic 
effects in the spectral functions presented in sections flV Bl 
and0 One obvious problem is that at high temperature 
both the number of data points and physical extent of 
the temporal direction becomes small, which make the 
reconstruction of the spectral functions difficult J13|. To 
some degree this problem can be overcome if one uses 
the information about the high energy behavior of the 
spectral function extracted at low temperatures, as was 
discussed in section IV. This issue will be discussed here 
in more detail. 

We start the discussion of the systematic effects for 
spectral functions corresponding to ground state, i.e. 
pseudoscalar and vector ones. In our method of analysis, 
the spectral function below plays a crucial role since 
we extract the high energy information from it, and we 
examine first the default model dependence of this spec- 
tral function. In Fig. E|we compare the pseudoscalar 
and vector spectral functions at 0.75Tc reconstructed us- 
ing the free massive lattice spectral function as a default 
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FIG. 12: The pseudoscalar (top) and vector (bottom) spec- 
tral functions at 0.75Tc reconstructed using the lattice and 
continuum default model (see text for further details) . 

Therefore, our way of using the information from the 
high uj part of the spectral function below Tc in the de- 
fault model for the analysis above T^ removes the arbi- 
trariness of the default model. To check this, in Fig. [T51 
we redo the analysis presented in Fig. 1101 where now the 
high UJ part was taken from the spectral function at 0.9 
Tc calculated with the free massive lattice default model 
with the corresponding mass. For comparison, we also 
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show the spectral function at 0.9 Tc recalculated with 
this default model. As one can see, the main features of 
Figs. El a) and b) are very similar to those of Fig. llOl al 
and b), respectively, so that one can essentially substi- 
tute them without having to change any of the discussion 
following Fig. E| 
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FIG. 13; The spectral functions for a) pseudoscalar and b) 
vector channels at different temperatures for set III, where 
the default model incorporates the high energy part of the 
spectral function at 0.9 Tc obtained using the free massive 
lattice spectral function (see text for further details). 

Another question to be addressed is to what extent 
the temperature modification of the spectral functions 
seen in Sec. |3 and Fig. El are real physical effects, as 
opposed to the offshoot of the inability of the MEM to 
correctly reproduce the spectral function from the small 
number of data points and small temporal extent avail- 
able at higher temperatures. To answer this question wc 
compare in Fig. ^Ithe spectral functions for set III at 
different temperatures constructed using the same num- 
ber of data points. In the three panels of Fig. El a) the 
pseudoscalar spectral functions at different temperatures 
are constructed with the same number of data points as 
are available for the highest temperature in the panel, 
omitting points from the center for the lower tempera- 
tures. N(jata here refers to the number of points after 
periodic folding. Though the figure shows that the de- 



tails of reconstructed spectral functions depend on the 
number of data points used in the analysis, it is clear 
that the fact that substantial modification of the spec- 
tral function takes place at 2.2hTc and the ground state 
is "dissolved" at 3Tc are not artifacts of the limited res- 
olution. Fig. El^^) also show very clearly the absence of 
any significant change in the pseudoscalar channel upto 
1.5 Tc- Similar comments also hold for the vector chan- 
nel in Fig. Elb). In this case, we do see the mass shift 
at 1.5 Tc discussed before irrespective of Njataj but the 
exact amount of the shift differs. 
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FIG. 14: The a) pseudoscalar and b) vector spectral func- 
tions at different temperatures reconstructed using the same 
number of data points at each temperature. 

In Fig. El we do a similar analysis for set II, where 
we reconstruct the spectral function at different temper- 
atures using the same physical extent as available for 
1.5 Tc- This figure shows that the seeming reduction of 
the peak strength at 1.5 Tc in Fig. |Hlwas probably due 
to the small number of data points available there, and 
the pseudoscalar and vector peaks are not substantially 
changed upto 1.5 Tc, as one would expect from Fig. Ela) 
and the results for set III. 

Now let us discuss the scalar and axial vector channels. 
In Fig. El we show the spectral function in the scalar 
and axial vector channels at O.OTc and set III using the 
free massive lattice spectral function as a default model. 
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FIG. 15: Spectral functions for set II for pseudoscalar and 
vector with Ndata = 6 data points. 



Also shown there are the spectral functions for free con- 
tinuum default model from section IV for comparison. 
As one can see the spectral functions show a stronger de- 
fault model dependence than the pseudoscalar and vector 
ones, though the main features of the spectral functions 
are independent of the default model. The stronger de- 
fault model dependence is probably due to the fact that 
these correlators are noisier than the pseudoscalar and 
vector ones. The default model dependence of the scalar 
and axial vector spectral functions was found to be even 
stronger for set II at 0.75Tc. If one uses the free massive 
spectral functions as default model at O.TSTc a large re- 
duction in the peak height of the Xc state can be seen, 
though the position of the peak is roughly correct. This 
problem is also due to limited statistics. It turns out that 
the probability P[(t|q;] = exp(L — aS*), which is expected 
to be strongly peaked around some amax ideally (for large 
statistics), has a long tail at large a values in this partic- 
ular case, so that the peak structure seen around amax is 
smoothened considerably when one averages over a. 

Finally, we also want to address the issue whether the 
serious system modification, possibly dissolution, of the 
Xc states at 1.12 Tc seen in Fig. [7|is an artefact of the 
smaller number of data points at this temperature. In 
Fig. El we show the spectral functions in scalar and ax- 
ial vector channels at 0.75Tc and 1.12Tc using Ndata = 8 
data points and the same temporal extent. The figure 
clearly shows that the large change in the spectral func- 
tion between 0.75Tc and 1.12Tc is a physical effect and 
not due to the small number of data points at the higher 
temperature. 

So far all data above Tc have been analyzed using the 
high energy part of the corresponding spectral functions 
below Tc- For set III at 1.5Tc we still have Ndata = 12 
data points and not too small temporal extent. Therefore 
here we can also reliably reconstruct the spectral func- 
tion without providing precise information for the high lu 
region. Figure ITSl show the spectral function for J/ip, re- 
constructed using both the free massless continuum and 
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FIG. 16: Spectral functions for scalar (top) and axial vector 
(bottom) channels reconstructed using the continuum and the 
lattice free spectral function as a default model. 



the free massive lattice spectral functions as the default 
model. Although some dependence on the default model 
is clearly seen, the strong ground state peak is pretty sta- 
ble and provides further evidence that the J/ip survives 
well into the deconfined gluon plasma. 



VII. SPATIAL CORRELATORS AND 
SCREENING MASSES 

Most of the early studies of the finite temperature 
meson spectrum concentrated on studies of correlators 
in one of the spatial direction (for convenience, we will 
use the z-direction in what follows). The corresponding 
masses are usually referred to as the screening masses 
[S^l- The spatial correlators are related to the same 
spectral function, but the relationship is a bit more com- 
plicated: the momentum-projected correlators in the z 
direction are given by [35l | 



G{iuJn,PT,z) = / ^ 



dpo 



2paa{pa,pT,Pz 



(24) 
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FIG. 17: Spectral functions from set II for scalar and axial 
vector with Ndata = 8 data points at different temperatures. 
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FIG. 18: The vector spectral function at l.STc and set III 
using the free continuum and lattice spectral function as a 
default model. 



where pt is the transverse momentum in the xy-plane 
and LOn is the Matsubara frequency. The screening mass 
is obtained from the correlators projected to zero trans- 
verse momentum and zero Matsubara frequency. 

While the screening masses, in general, are not directly 
related to the spectrum of the finite temperature sys- 
tem, it is possible to extract important informations from 
them, in particular for examining particular model spec- 
tral functions. In the presence of a stable bound state, 
that gives a contribution ^ S{pq — — m^) to the spec- 
tral function, it can be easily seen from Eq. H24|l that the 
screening mass is identical to the pole mass of the state. 

Numerically, the study of screening masses is much 
simpler because of the large extent available in the z- 
direction, and one can obtain the screening masses from 
standard one-exponential fits to the long distance part of 
the correlator. Fig. 1191 summarizes the temperature de- 
pendence of the vector and pseudoscalar screening masses 
for the different sets studied by us. We have calculated 
them from one-exponential and two-exponential fits, and 



in all cases the two agree. In order to show the rela- 
tive changes in screening mass with temperature we have 
normalized the masses with the zero-temperature masses 
from Tabic IIVI which, as mentioned in Sec. IIIII were 
obtained from the screening masses from the lowest tem- 
perature available to us for each set. So the part below 
Tc in Fig. 1191 is just the normalization except in Set I. 
From the figure, one makes these observations: 

1) The screening mass does not change upto Tc (or at 
least, till 0.93 Tc) for either channel. This observation 
is in agreement with what is seen for light quarks and 
glueballs. We also see that the screening masses below Tc 
are in agreement with the pole mass extracted from the 
temporal masses, wherever such an extraction is possible 
reliably. 
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FIG. 19: Screening masses for the a) vector and b) pseu- 
doscalar channels for the different sets, normalized by the 
zero-temperature mass (see text). 

2) The screening masses show a definite temperature 
effect above Tc, increasing steadily with temperature for 
both the channels. This effect is in sharp contrast to the 
temporal correlators, which indicated little temperature 
effect upto « 1.5Tc, in particular for the pseudoscalar 
channel. Since for the pseudoscalar channel, we do not 
see any shift in the ground state peak position certainly 
upto 1.5 Tc, and possibly upto 2.25 Tc, Fig. [Hb) in- 
dicates that the pole and screening masses differ from 
each other, even though a bound state is present in the 
spectral function. For the vector case also, even though 
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a small shift in the pole position may be present at 1.5 
Tc, the shift of screening mass in Fig. E| a) is much 
larger and indicate a similar difference between the pole 
and screening masses. One possibility that one has to 
consider is that while there is a bound state in the spec- 
tral function above Tc, it is not a S function state. We 
investigated the possibility of the screening mass change 
being caused by the state gaining a width above by 
using a Breit-Wigner form as found in Ref. ^3 j but our 
conclusion is that the change in screening mass is too 
large to be explained this way. 

Another possible way that a screening mass can change 
differently from the pole mass, even in the presence of a 
bound state, is if the dispersion relation changes in the 
thermal medium "s^. Since in a thermal medium, the 
Lorentz symmetry is lost, in general one would expect an 
energy-momentum asymmetric self energy term, leading 
to a dispersion relation 

Lo^ip, T) = m{Tf +f + n(p, T). (25) 

This will, in general, alter the screening mass from the 
pole mass. For small momenta, one can make the sim- 
plifying assumption that the modified dispersion relation 
can be approximated by a temperature-dependent mass 
and a temperature-dependent "speed of light" [s^ 

u;\p, T) = m'iT) + A^{T)f (26) 

For a delta-function state in the dispersion relation, Eq. 
gives 

(T) = TOpoic('r)/A(T). (27) 

Such a modified dispersion relation can explain the dis- 
crepancy between the pole and screening mass if, for ex- 
ample, A{T) < 1 at 1.5 Tc. 

This scenario can be checked by constructing the spec- 
tral function from the temporal correlators for non-zero 
momenta. Eq. (|26|l with A{T) < 1 will lead to a less- 
than-relativistic shift in the spectral function peak with 
spatial momenta. Figure |^ a) shows the ground state 
peak for the vector channel below and above Tc for set 
III, and Fig. [201 b) shows the same for the pseudoscalar 
channel. The vertical bands in both the figures show 
the expected peak positions for the non-zero momenta, 
assuming a relativistic shift from the peak at zero mo- 
menta. This was obtained by using the lattice form of 
Eq. (Uni with A{T) = 1, 

2 

coshcj(p) - 1 = ^ (1 - cosp,) + ^, (28) 

1=1,2 

where all the variables uj,pi and m are in units of inverse 
lattice spacing. The width of the band reflects the un- 
certainty in the zero momentum peak position due to the 
finite bin size. The shift of the peak position at 0.9 is 
seen to be consistent with the expected relativistic shift 
in both the channels, in agreement with previous results 



that the Lorentz symmetry is approximately valid upto 
close to Tc Hmf. At 1.5 T^, one clearly sees a much 
smaller shift in the peak position, confirming the medium 
modification of the dispersion relation discussed above. 
A more quantitative analysis, to check whether the shift 
in the screening mass is completely due to such a mod- 
ification of the dispersion relation, cannot, however, be 
conclusively made due to the rather large uncertainty in 
determining the peak position from a maximum entropy 
analysis. 
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FIG. 20: Spectral functions for the a) pseudoscalar and b) 
vector channel for a few non-zero momenta, below and above 
Tc. The momenta are given in units of j^, where La is the 
transverse dimension of the lattice. The vertical bands show 
the expected peak positions for the non-zero momenta, as- 
suming a relativistic dispersion relation. 



VIII. SUMMARY AND DISCUSSION 

In this work, we have conducted a detailed lattice study 
of the properties of the IS {r]c and J/ip) and IP {Xc 
and Xc) charmonia in hot gluonic medium, by studying 
the thermal correlators in Euclidean time and extract- 
ing spectral functions from them. We show that one 
can extract useful information about medium modifica- 
tions of charmonia above Tc by comparing the correla- 
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tors with those reconstructed from the spectral function 
below Tj,. We also show that by providing correct infor- 
mation about the structure of the high energy part of the 
spectral function, one can reliably extract the low energy 
part of the spectral function at high temperatures from 
correlators measured at a limited number of points. 

In the hadronic phase, we find that the properties of all 
the states remain unchanged at least upto temperatures 
of 0.93 Tc as one approaches the deconfinement transition 
temperature from below. This is similar to quenched 
studies of light mesons and glueballs, and may not be 
true for full QCD. 

As the temperature of the medium crosses Tc, signifi- 
cant modifications of the x states are observed. The tem- 
poral correlators show large enhancement, and already 
at 1.1 Tc, the ground state peaks are not observed in 
the spectral function, indicating that these states may 
have dissociated already at this temperature. The IS 
states behave quite differently. Little change is seen in 
the correlators as one crosses Tc, and the spectral func- 
tion shows a significant ground state peak till quite high 
temperatures. No significant reduction of either the peak 
strength or the mass of the state is seen at least upto 1.5 
Tc. At higher temperatures, the peak weakens, though 
at 2.25 Tc one still finds a significant peak. (The appar- 
ent contradiction of this result with Ref. 'l4| has been 
discussed in Sec. 0) Finally, at 3 Tc we do not ob- 
serve any significant peak any more, indicating that the 
state has broadened and weakened so much that it is not 
meaningful to treat it as a resonance any more. 

These results have direct phenomenological implica- 
tions, as p{u)) for the vector current is connected to 
the dilepton rate. Let us consider a hypothetical heavy 
ion collision experiment which forms equilibriated gluon 
plasma, and whose temperature can be increased gradu- 
ally. As one crosses Tc, one expects to see a reduction in 
the J /tp peak in the dilepton channel caused by the disap- 
pearance of the excited states. (We would like to remind 
here that about 40% of all J/ -0 produced in hadron colli- 
sions are indirect and come from Xc and ip' states.) The 
remaining J /ip peak then stays stable as one increases 
the temperature further, showing little further reduction 
at least upto 1.5 Tc. On further increasing the temper- 



ature, the peak starts to "melt" gradually by weakening 
and possibly broadening. At 2.25 Tc we see ~ 25% fur- 
ther reduction in the integrated strength of the J/ij) peak 
(for our quark mass). This process continues as one in- 
creases the temperaure further, and by 3 Tc the J /ip peak 
is indistinguishable from the background. 

It is reasonable to expect that the qualitative picture 
remains the same as one introduces dynamical quarks. 
Below Tc, there may be changes in the meson proper- 
ties due to the activation of the thermal pions, and it is 
possible that the x states dissociate even before Tc. For 
J /tp itself, such thermal pions will have little effect, and 
above Tc, the gradual melting picture described above 
is likely to be still true qualitatively, with the thermal 
quarks causing only small quantitative changes. 

These results are quite different from the earlier po- 
tential model studies, which predicted the IS charmonia 
to be dissolved at ~ l.lTc Since the appearance of 

our earlier studies fl^ some of the potential model stud- 
ies have been reanalyzed. It was pointed out in Ref. [3^ 
that Ref. 2] does not take into account a possible rise in 
the strong coupling constant near Tc which will lead to a 
bound J/ tp deep in the plasma. Also a very thorough and 
detailed study of the color singlet potential conducted in 
Ref. shows that the color singlet free energy and 

potential above Tc (upto a few times Tc) behave very 
differently from what was anticipated in Ref. |^ on per- 
turbative grounds jJ9j, and will probably give rise to a 
considerably higher dissolution temperature for the J/ip- 
Be that as it may, it is not clear that a potential model 
study can capture all the physics going into the medium 
modification of J/tp like the collision broadening due to 
thermal giuons. It will certainly be interesting to have 
an analysis based on potential models or other similar 
studies that can reproduce all the features of the direct 
studies by us (and others), which also rules out any sub- 
stantial reduction of the mass of the J/4' and rjc above 
Tc from its zero temperature mass. 
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